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Abstract. We prove a GUE central limit theorem for random variables with finite fourth moment. We 
apply this theorem to prove that the directed first and last passage percolation problems in thin rectangles 
exhibit universal fluctuations given by the Tracy- Widom law. 



1. Introduction 

In the last few years, it has become clear that random matrix theory is intimately related to a variety of 
questions arising in physics, statistics, combinatorics, representation theory, number theory, and probability 
theory (see e.g. |12ll51 [Tnill9p . It is the fluctuations of random matrix ensembles which are common to the 
specific problems from each of these areas. In the hope of understanding the universal nature of the random 
matrix distributions, it is natural to search for central limit theorems for which these distributions are the 
limiting objects. 

The following is a Gaussian Unitary Ensemble (GUE) central limit theorem. 

Theorem 1 (GUECLT). Suppose that {Xl}°°^^ is a family of independent identically distributed random 
variables such that EX/ = 0, E|X/P = 1, and ElxU^ < oo. Let 



(1) L{N,k)= sup J2{ E ^'}' 

k ij 

(2) R{N,k)= inf y{ y\ ^'}- 

Q=io<h<...<ik=N ^ ^ ^' 
j=l l=lj-l+l 

If k,N ^ oo such that k = o{N°'), a < jg, then 

(3) (^^2Vk)ki=^FauE, 



(4) (-^-2Vk)ki=^FauE 



R(N,k) 

— Z V K I K o =? 

where Fqjje is the GUE Tracy-Widom distribution. If {Xf}°°^-^ are independent identically distributed 
Gaussian random variables satisfying the above conditions, then (PJ and (@J hold if k = o{N"), a < |. 

The GUE Tracy-Widom distribution function ^S] is given by 
(5) Fgue{x) = exp< - (s - x)q'^{s)ds }, 



where q{x) solves the Painleve II equation, 

(6) g" = 2q^ + xq 

subject to the condition that q{x) ^ Ai(x) as x +00; Ai(x) denotes the Airy function. The function 
Fque is the limiting distribution function for the largest eigenvalue of the Gaussian Unitary Ensemble as 
the dimension of the matrices grows to infinity (see H18|l 'l. 



Theorem is intimately related to the directed first and last passage percolation problems. Consider 
the N X N lattice and a set of associated independent identically distributed random variables {^f 
satisfying EXf = /i, E|X/p - = cr^, and E|Xf I"' < oo. An up/right path tt from the site (1, 1) to the 
site {N,k) is a collection of sites {{h, jk)}h^i~^ satisfying = (1,1), {iN+k-i, iN+k-i) = {N,k) and 

(ifc+i, jfc+i) — {ik,jk) is either (1,0) or (0, 1). Let (1, 1) y {N, k) denote the set of such up/right paths. The 
directed first and last passage times to [N, fc) S N x N, denoted by [N, k) and L'-{N, k), respectively, are 
defined as 

(7) Lf{N,k) = min V Xl, 

Te(l,l)/-(Ar,fc) 

(8) LHN.k) = max V X/. 

^e(i,i)/-(jv,fc) 

If xl is interpreted as the time to pass the site («, j), L^{N, k) and L^{N, k) represent the minimal and the 
maximal time to travel from the site (1, 1) to [N, k) along an admissible path. Since the directed last passage 
percolation time can be viewed as the departure time in queuing theory (see e.g. [7]), the result below also 
applies to queuing theory. In addition, Corollarv ll.il also applies to the flux of particles at a given site in 
the totally asymmetric simple exclusion process (see e.g. 221)- Theorem Q] implies the following. 

Corollary 1.1. Suppose that {Xf^f^^^ is a family of independent identically distributed random variables 
such that EXl = ^, E|X/P - = cr^, and E\Xi\^ < oo. For any seM., 

^/LUN,k)- n(N + k-l)-2aVNk \ 
9 lim P ' ' — ^ ■ H <s] = FauEis), 

(10) Jim ak-^/^Nh "V ^ ^-Fgue{-s), 



where k = o{N°') and a < jg- If {Xl}°°^^ are independent identically distributed Gaussian random variables 
satisfying the above conditions, then ^ and H10(l hold when k = o{N"), a < |. 



The restriction k = o{N"), a < in the above corollary seems to be a technical matter. It is believed 
that such a central limit theorem (after possible changes in scaling) should hold as iV, fc ^ oo with no 
restrictions on their relative rate of growth. When N and k are of the same order, such a limit law for was 
obtained for geometric and exponential random variables by Johansson using Young tableaux theory in 
combinatorics and techniques from random matrix theory. However, analyzing the last passage percolation 
time for general random variables in arbitrary scaling regimes is still an open question. 

Even the determination of the 'time constants' as iV, fc — > oo for general random variables is a challenge. 
Glynn and Whitt [7j proved that for k = o(7V), 

(11) lim m^)^^^^ 



P-almost surely for some a, where a was proved to be equal to 2a by Seppalainen (see also JJ). Note 
that the definition of the last passage percolation, L''{N,k), is symmetric in N and k. Hence, (|ll|l can not 
be true when A'' and k are of the same order. For exponential random variables of mean 1, an interpretation 
of a theorem of Rost ^B] yields that 

(12) hm ^'^[""]'[^"]) =(Vi+Vy)^ 

n — >oo fi 



For slightly different choices of admissible paths, a few more explicit random variables are shown to have the same type of 
limit law (see e.g. ITIISI 'I 
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Johansson |5] (see equation (1.4)) proved that for geometric random variables of parameter q, 

L'{[xn],[yn]) q{x + y) + 



(13) hm 

n^oo n 1^9 

These two cases seem to be the only explicitly computed examples when N and k are of the same order. 
Note that the directed first and last passage times can also be written as 

(14) Lf{N,k) = inf H XI 

(15) L'(7V,fc) = sup ^( Xi 

When k is fixed and N ^ oo, the last passage time is obtained along a path that lies on the first level, j = 1, 
for a certain 'time', [tiiVj , then jumps to the second level, j = 2, and stays there until the 'time' [t2AfJ , and 
so on such that the total sum is maximized. As the sum of random variables on each level is 'asymptotically 
equal' to a Brownian motion after proper centering and scaling, by the Donskcr's theorem one can imagine 
that 

(16) ^L(^Xl^_if^ ^ ^,.(1) := 3up ^(i?,(i,)-5,(t,_0), 

where Bj(t), j ~ l,...,fc, are independent Brownian motions. This was proved in [7| in the context of 
queuing theory. 

The distribution function of 1)^(1 ) is difficult to compute since Dk{l) is a complicated and unusual 
functional of k Brownian motions. However, as mentioned earlier, for some special random variables. Young 
tableaux theory in combinatorics can be applied to explicitly compute the distribution of LK By studying a 
special random variable, Baryshnikov [3j and Gravner, Tracy and Widom |H] (see also |14II13| ) showed that 

k 

(17) p(^fe(i)<.)=^ r ... r n i^^-o^ne-^^^^rfe,, 

where Zk = l!2! . . . fc!(27r)''^/^. As the integrand on the right-hand-side is the density of the eigenvalues 
^1, . . . , ^fc of fc X fc GUE random matrix, this identity implies that -Dfc(l) has precisely the same distribution 
as the largest eigenvalue, Cmax(fc), of fc x fc GUE random matrix. It is well-known in random matrix theory 
that (see e.g. piT^) 

(18) lim P((Cmax(fc) - 2Vfc)fci/^ <s)= FauEis). 



Therefore, 



(19) lim lim f( | ^'(^^fc)j/^^ „ 2Vk]k'/' < s] = Fgue{s). 

The content of CoroUarv ll.ll is that the two limits can be taken simultaneously as long as /c = o(iV"), a < ^i- 
Note that in the centering fj,{N + k'~ 1) of 0, the term fJ,{k — 1) accounts for the 'up movements' of the path 
while fiN comes from the 'right movements' of the path. The term fi(k — 1) is negligible when k ~ o(N~), 
but it is retained in the formula in order to emphasize the symmetry of N,k. 

Theorem n is a central limit theorem whose limit is the Tracy- Widom GUE top and bottom eigenvalue 
distributions. The joint distribution of the top and bottom n eigenvalues is also a natural limit of a central 
limit theorem which depends on the O'ConncU and Yor representation of the /c-dimensional Dyson process |14l 
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IB] . Following [Hini, let Co([0, 1],R) be the space of continuous functions / : [0, 1] ^ R with /(O) = 0. 
For f,g£ Co([0, 1],R), define 

(20) / ® g{t) = inf ;/Gs) + 9{t) - g{s)], 

(21) /0 5(t)= sup [/(s)+.g(t)-.g(s)]. 

0<s<t 

The order of operations is from left to right. Define a sequence of mappings Ffc : Cq ( [0, 1] , R)'' Co ( [0, 1] , R)*^ 

by 

(22) r2(/,5) = (/®5,5 0/), 
and for k > 2, 

(23) rfe(/i,. ..,/,) = (/i®/20/3 0---®/fc, 

rfe-l(/2 /l, /3 (/l /2), /fe (/l • • • fk-l))). 

O'Connell and Yor proved that if Bi, are independent standard one-dimensional Brownian motions, 

Ffc(i?i, Bk) has the same distribution on path space, Co([0, 1], M)''", as k one-dimensional Brownian motions 
starting from the origin, conditioned (in the sense of Doob) never to collide. It is well known that this process 
can be interpreted as the fc-dimcnsional GUE Dyson eigenvalue process. The first coordinate of F^ is the 
smallest eigenvalue, the second coordinate of F^ is the second smallest eigenvalue, and so on. The proof of 
Theorem ^ immediately implies the following theorem once one notes that the first n coordinates of F^ are 
Lipshitz functions with a fixed Lipshitz constant independent of k. We omit the details. 

Theorem 2. Suppose that {X^^}f°^^ is a family of independent identically distributed random variables 
such that EXil = 0, E|X./|2 = 1, and < oo. Let s%{i) = -^J2l=i^i- As N,k oo such that 

k^O{m),a<f-^, 

(24) (F,( )i(A^),...,r,( UN))=^FSuE, 

where Fqjj^ denotes the limiting joint probability distribution of the bottom n eigenvalues of the GUE. By 
symmetry, a similar statement holds for the top n eigenvalues. 

Proving a GUE central limit theorem and analyzing certain scaling regimes for the directed first and 
last passage percolation problems is the main focus of this paper. It is interesting to seek central limit 
type theorems for the Gaussian Orthogonal Ensemble (GOE) and the Gaussian Symplectic Ensemble (GSE) 
Tracy- Widom distributions. The arguments of this paper can be adapted to prove such theorems. The 
authors will address this topic in several forthcoming papers. The authors will also describe the fluctuations 
of restricted Plancharel measures from this point of view. 

Note. While completing this paper, the authors learned that T. Bodincau and J. Martin had announced a 
similar result 0] . Both 0] and this paper use the idea of approximating random walks by Brownian motion 
in order to analyze the last passage time. As opposed to the authors use the Skorohod embedding 
theorem to couple random walks to Brownian motion. Bodineau and Martin use the Komlos, Major and 
Tusnady approximation theorem (which produces tighter bounds) to achieve such a coupling. The authors 
arc grateful to Rongfcng Sun for bringing (5 to their attention. 

Acknowledgments. The authors thank S.R.S. Varadhan for his crucial comments at an early stage of this 
work. The authors thank G. Ben Arous, P. Deift, and R. Sun for useful discussions. The authors also thank 
J. Martin for constructive comments on the first version of this paper. The work of Baik was supported in 
part by NSF Grant #DMS-0350729 and the AMS Centennial Fellowship. The work of Suidan was supported 
in part by a NSF postdoctoral fellowship. 
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2. Proofs 



Theorem n and Corollarv ll.ll arc proven in this section. 

Proof of Theorem^ The proof is based on embedding general random walks in Brownian motions by means 
of the classical Skorohod embedding theorem 

Theorem 3 (Skorohod Embedding Theorem). Let {-Bt}tgK+ be a one- dimensional standard Brownian mo- 
tion and X be a real valued random variable satisfying EX = 0, EX^ = 1, and EX* < oo. Then, there exists 
a Brownian stopping time, t, such that Br is distributed as X , Et = EX^ = 1, and Er^ < 4EX'*. 

Let {Xi}°^i be a sequence of independent identically distributed random variables satisfying the assump- 
tions of Theorem|2| There exists a sequence of independent identically distributed positive random variables 
{Ti}°^^ such that Bri+...+T„ is distributed as Xi + ... -\- Xn and {Bri+...+Tr,+i — ^Ti+...+r„}$^Lo is a sequence 
of independent identically distributed random variables with distribution Xi. This follows from Theorem 
and the strong Markov property of Brownian motion. 

To prove uniform estimates in time for the difference of embedded random walks and the scaled Brownian 
motion in which the walks are embedded, the following two real valued processes are useful: 



(25) SN{t) 

(26) BN{t) 



Sn+...+rL«,j + {Nt- [Nt\){Br,,,, -Br,,,,) 



N 



B 



Nt 



Let C([0, IJjM'^) be the space of continuous fc-vector valued functions on the unit interval equipped with the 
sup norm, || • ||fc: for / G C([0, 1],M''), ||/||fe = sup^gjo^ij \fj{t)\. Theorem[T]is proven in two steps. 

The first step involves establishing probabilistic estimates on differences of the above processes. The second 
step involves a Lipshitz bound on the relevant Brownian concatenation operations and the application of the 
quantitative estimates established in the first step. 

Two estimates will be used in the proof of TheoremQ] Let Yi = Ti — 1 and Z^^ = J2n=i ^n- Note that 
Zf and [Z^Y are submartingales. By the Doob martingale inequalities [TB). 



»f sup |Zf I > /3n) 



N\2 



< 



(27) < 



nEJ=l^^? _ Erf- (En )^ 
Var{Ti) 



For Pn = iSt where c > and A < i this inequality implies 



(28) Pf sup " '- > — T ) < 

This is the first estimate. 

The second estimate concerns the modulus of continuity for Brownian motion. If p < ^, then for some 
appropriately chosen constants A,i^>0, 

(29) P( sup \BN{t) - Bn{s)\ > ^) < Ae-''^ 

He[0,l], \t-s\<cN-^ ^^^^ 

uniformly in A^, c > 0. This estimate is a consequence of Brownian scaling and standard estimates for the 
maximum of Brownian motion. The first part of the proof of Theorem Q is complete. 
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Consider the functions Gk ■ C([0, 1],R''") M and : C([0, 1],R''") ^ R defined by 



k 

(30) G,(/) = inf - /.fe-i))' 



k 

O=to<ti<...<t& = l ■ 



(31) Gfc(/) = _ ^ sup ^J2yAh)~fAh~i))- 



Lemma 2.1. iJoi/i Gfc and G^ are Lipshitz with Lipshitz constant 2: 

(32) |Gfc(/)-Gfe(5)| <2||/-g|U and |Gfe(/) - Gfe(g)| < 2||/ - g|U, 
/or aZ/ e G([0, f], R'^). 

proof: Let f,g& G([0, 1], R*-') satisfy = e, \\fi-gi\\i = e^, andX;*Li = e where / = {fi, fk), g = 

{gi,-9k)- _|Gfc(/) - Gfe(g)| < X;J=i \Gk{.fi,-;fi-i,9i,-;9k) - Gk{fi, fi, gr+i, ■■■gk)\_hy the triangle in- 
equality. |Gfe(/i,...,/i_i,5i,...,5fe) - Gfe(/i, ...,/,, 5i+i,...gfe)| < 2e,;. Thus, \Gk{f) - Gk{g)\ < 2j|/ - 
Since the same argument holds for |Gj.(/) — Gi.{g)\, the proof of the lemma is complete. O 

By Skorohod embedding, for each A'^ > construct a sequence of pairs of random walks embedded in 
standard Brownian motions, {(5';^(i), as in H25() and H26() . Lemma r2.1l and the estimates on the 

differences of the Skorohod coupling imply that if A^, fc ^ oo and k ~ o{N°'), a < jg, then 

(33) k^Gk{Sl„...,S'^)~Gk{Bl,,...,B%)) ^0, 

(34) k-e (G,(5^, S%) - GkiBl,, B%)) ^ 0, 
in distribution. More explicitly, if fc = [A^"J , then 

_ _ 

F{k^Gk{S'^,...,S%)-Gk{Bj„...,B%))>e) < p( ^ - B^Hi > 



1=1 



2N6 



(35) < ^E||5i,-Bi.||i. 



Note that 



E||5^-i3^||i = m\Sl, - Bl,\\, > s)ds ^ — p(^\\s},-Bl,\\,>—yu 

= — / P(||5i.-Bi.||i > — ; sup \Y,r,~A>uN'-^)du 
^vjo ^ i<i<N ~; ^ 



(36) +— / P ||5i.-S^||i>— ; sup |5Jr,-z|<«7Vi-Mdw. 

JV^ Jo ^ l<i<Ar ^ 

If A < i, then H28|) implies that the integrand of the first integral on the right hand side of (|36|l has an 
integrable tail. If the {Xil}°°^^ are Gaussian, this term is identically since ti = 1 P-almost surely; in this 
case, there is no constraint on A. The integrand of the second integral of H36() is dominated by 

(37) Pf sup \BN{t) ~ Bn{s)\ > 



te[0,l],\t-s\<uN->' N''^ 

(|29|l implies that if p < -j, the second integrand has exponentially decaying tails uniformly in N . A necessary 
condition for the right hand side of H35|l to vanish as ^ oo is a < fp. If {Xl}f^^^^ are Gaussian, then 
p < i is sufficient to make the second integral of (I36|l finite; thus, if a < |, the right hand side of (|35|) 
vanishes as iV ^ oo. If {Xl^fj^^ are not Gaussian, then a < |p and 2p < A < i imply that a < ^ is 
sufficient to make H36(l vanish as iV oo. 
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As mentioned in the introduction, the theorems of Baryshnikov [2 and Gravener, Tracy, and Widom [5] 
state that Gfc(i?]^, is distributed as the largest eigenvalue of the GUE of fc x fc matrices. By (I18|) 

and 



(38) 



F, 



GUE, 



in distribution. By H34|) and the symmetry with respect to multiplication by —1 of Brownian motion and 
the Dyson process 



(39) 



k-^{'Gk{S],,...,S%) ^2Vk) Fgue. 



Since GkiSj^, S'^) and Gj.(S']^, S'^) are distributed as L{N, k) and R{N^ fc), respectively, the proof of 
Theorem n is complete. □ 



It is clear from the proof of Theorem^that when fc >> o(iV?) the 'up movements' of the paths contribute 
to L^{N, fc). Hence, one needs to take those terms into account in order to prove a central limit theorem for 
general fc scaling. The arguments in this paper do not seem to extend to the general case. 

The proof of CoroUarv 1 1 . II remains . 

proof of Corollary With no loss of generality, assume that EX/ = and E|X/p = 1. The definitions of 
L'-{N, fc) and L{N, fc) imply the following deterministic estimate: 



fe-i 



(40) fcs inf i^yX^+^]<k-o 

l<n<--<ik-l<N ^ ^ J 

Note that 



PiN,k)-L{N,k) ^^^1 



k-l 

"(ki sup -^yxi+' 



< 



N 



ke 



sup I 1 



fc-1 



sup 



ke 



< 



< 



ki 
Wn JO 



(41) 



ke 

ey/N JO 

ke 

eNi Jo 



E( sup \Xl\) 

{l-P(|X]^| < s)^}ds 

"(iXll > sN-j ids 



1 - 



1 - 



Now split the last integral into two pieces; one over the interval [0, 1] and the other over (!,{»). Using the 
basic estimate < P^|Xj| > sN^^ < 1 for the first integral and the Chebyshev inequality for the second 
integral, H41|l is less than or equal to 



(42) 



eNi Jo 



— n 



1 



1 



E\Xl\ 



yds < 



ke 



fce 



ds. 



If iV, fc ^ oo such that fc = o{N"), a < jg, then the right hand side of 1)42(1 tends to 0. The inf can be 
treated similarly. Theorem^and H41() . (|42|l imply CoroUarv 1 1.1 1 for the non-Gaussian case. If {XfJ'^j^^ are 
Gaussian, then the arguments from the proof of Theorem ^ apply. □ 
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